Titanate-based high-entropy perovskite oxides relaxor ferroelectrics

Different combinations of monovalent and trivalent A-cations in high-entropy perovskite oxides (HEPOs) were investigated. The multicomponent (A′0.2A″0.2Ba0.2Sr0.2Ca0.2)TiO3 (A′ = Na+, K+, A″ = Bi3+, La3+) perovskite compounds were successfully synthesized by solid-state reaction method persisting average cubic perovskite phase. The trivalent cation exhibited distinct effects on local structure, dielectric properties and relaxor ferroelectric behavior. Highly dense ceramics (> 95%), high dielectric constant (~ 3000), low dielectric loss (~ 0.1), and relaxor ferroelectric characteristics were obtained in the compound containing Bi3+. The La3+ containing compounds revealed lower dielectric constant, higher dielectric loss and linear dielectric behavior. The effect of monovalent cation on the dielectric properties was minimal. However, it affected relaxor ferroelectric behavior at elevated temperatures and conduction behavior at high temperatures. The (K0.2Bi0.2Ba0.2Sr0.2Ca0.2)TiO3 ceramic maintained the relaxor ferroelectric behavior with low PREM at high temperatures suggesting more stable relaxor ferroelectric characteristics than that of the (Na0.2Bi0.2Ba0.2Sr0.2Ca0.2)TiO3. Moreover, between these two compounds, the homogeneous electrical characteristics could be obtained from the compound consisting of K + and Bi + at A-site. This study suggests that tuning the chemical composition, particularly choosing appropriate combination of mono/trivalent cations in high entropy perovskite oxides, could be the effective approach to develop high-performance relaxor ferroelectrics with the desired properties.


Experimental
(A′ 0.2 A″ 0.2 Ba 0.2 Sr 0.2 Ca 0.2 )TiO 3 (A′ = Na + , K + , A″ = Bi 3+ , La 3+ ) ceramics was prepared by a solid-state reaction method.Analytical grade reagents of Na 2 CO 3 , K 2 CO 3 , Bi 2 O 3 , La(OH) 3 , BaCO 3 , CaCO 3 , SrCO 3 , and TiO 2 were used as the raw materials.All starting powders were dried at 100 °C for 24 h to remove moisture.The starting powders were stoichiometrically weighed and mixed by ball milling method in ethanol by using zirconia grinding balls for 24 h.The slurry was dried by vacuum evaporation process with the speed of 90 rpm at 90 °C.The dried powders were calcined in air atmosphere at 975 °C for 2 h with a heating/cooling rate of 5 °C/min in a covered alumina crucible.The 2 wt% polyvinyl alcohol solution was added to the calcined powders by mixing using mortar and pestle until forming dried powders.They were then uniaxially pressed in a cylindrical die under 80 MPa for 2 min to obtain disc pellets or green pellets with 10 mm in diameter and about 1 mm in thickness.The ceramics were sintered at 1250 °C for 2 h with a heating/cooling rate of 2 °C/min in covered alumina crucibles.Finally, the sintered ceramics were polished to less than 0.8 mm in thickness and coated on both sides with a silver paste as electrodes.The coated ceramics were fired at 600 °C for 1 h with the heating/cooling rate of 5 °C/ min to eliminate organic content in the silver paste prior to electrical properties measurements.
Rietveld refinement method was carried out using GSAS-II software 12 to study phase purity and crystallographic structure properties of the ceramics by using of synchrotron X-ray diffraction data collected using Debye-Scherrer geometry with a strip detector (Mythen6K 450 (DECTRIS®) at BL1.1W MXT of the Synchrotron Light Research Institute (SLRI), Thailand.The samples prepared by grinding the ceramic samples into powder and annealed at 400 °C for 1 h were used for synchrotron XRD characterization.Raman spectroscopy measurements (XploRa PLUS, HORIBA) were performed by using a laser wavelength of 532 nm to investigate a short-range structure of the ceramics.X-Ray Fluorescence Spectrometer ((XRF), BRUKER S8 TIGER) was used to determine elemental composition of the ceramics.Shrinkage measurement was performed by measuring linear change in diameter between green pellet and sintered ceramic.Bulk density was determined with Archimedes' method.Scanning electron microscopy (SEM) were used to investigate surface microstructure of the ceramics (SEM, JEOL JSM-6480LV).Grain size of the ceramics was measured using the linear intercept method.For electrical characterizations, dielectric properties as a function of temperature (25-350 °C) were determined using an LCR-meter (HP model 4192A) at frequencies ranging from 10 to 1000 kHz.Ferroelectric system based on Radiant Precision High Voltage Interface was used to measure the polarization-electric field (P-E) hysteresis loops between 30 °C (Room Temperature, RT) and 150 °C by applying an electric field in the range of 10-100 kV/cm and a frequency of 1 Hz.Similarly, by using the same data from ferroelectric properties measurement, W REC , W LOSS , and η values were also calculated.The impedance analyzer (Solatron, 1260A) was used to test the complex impedance and complex modulus under frequencies from 0.1 Hz to 1 MHz, and a temperature range from 400 to 500 °C.

Rietveld refinement of synchrotron XRD and Raman Spectroscopy
Rietveld refinement analysis of synchrotron XRD patterns of NB, KB, NL, and KL ceramics were depicted in Fig. 1a-d.The data including R w (weighted-profile factor), χ 2 , GOF (goodness of fit), crystal structure, lattice parameter, and unit cell volume obtained from refinement were listed in Table 1.The R w , χ 2 , and GOF were calculated to evaluate the fitting quality of the experimental data in order to achieve the best possible fit to the experimental diffraction data.All ceramics showed small values of R w , χ 2 , and GOF, indicating that refined results agreed with cubic structure (Pm-3 m space group).This suggested that NB, KB, NL, and KL ceramics form an average structure of cubic-perovskite phase.The formation of single-phase perovskite structure with cubic structure agreed well with the calculated Goldschmidt tolerance (t), as shown in Eq. ( 4) 13 where RA was defined as the average ionic radii of cation at A-site, including the following: Na + (1.39 Å), K + (1.64 Å), Bi 3+ (1.36 Å), La 3+ (1.36 Å), Ba 2+ (1.61 Å), Ca 2+ (1.34 Å), and Sr 2+ (1.44 Å).R B was defined as the ionic radii of cation at B-site, which is Ti 4+ (0.61 Å).R O was described as the radii of O 2− ion (1.40Å) [14][15][16] .The calcu- lated tolerance factor for NB, KB, NL, and KL ceramics were 0.995, 1.012, 0.995, and 1.012, respectively.The KB ceramic exhibits larger lattice parameter and unit cell volume than that of NB, NL, and KL ceramics, respectively.It indicated that the different mono-/tri-valent cations at A-site lattice of ceramics affect the structural formation due to ionic size difference of the A-cations.
The local structure of NB, KB, NL, and KL ceramics was determined by using Raman spectroscopy by identifying the vibration modes that persisted in the oxide structures.As seen in Fig. 1e, the four distinct bands corresponded to the main Raman modes of perovskite ceramics including 100-200 cm −1 , 200-450 cm −1 , 450-700 cm −1 , and above 700 cm −1 were observed for all compositions.The Raman spectrum of the prepared (Bi 0.5 Na 0.5 )TiO 3 ceramic was included as reference.The first band between 100 and 200 cm −1 correlated to the vibration modes of A-site cations (i.e., Na/K-O, Bi/La-O, Ba-O, Ca-O, and Sr-O bonds) in the perovskite structure.The second band in the range of 200-450 cm −1 correlated with B-site vibrations (Ti-O bonds).The third band of 450-700 cm −1 related to the stretching vibrations of TiO 6 octahedral.The A and E overlapping bands in the final region above 700 cm −1 related with oxygen displacements.The NB and KB spectra displayed distinct characteristics between 100 and 200 cm −1 than that of the NL and KL spectra.It could be attributed to different cation disorder at A-site caused by Bi and La ion resulting in different bonding characteristics of A-O bonds.Moreover, according to literature, the presence of hybridization of 6s 2 bismuth lone pair and oxygen p-orbital played a vital part in local structure of bismuth sodium titanate perovskites 17 .The weaker A-O bond could be induced by introducing La onto A-site and by increasing temperature to depolarizing temperature, suggesting loss of the hybridization of 6s 2 bismuth lone pair and oxygen p-orbital [17][18][19] .A clear observation in distinct A-O bands between NB/KB and NL/KL suggested that the presence or the lack of hybridization of 6s 2 bismuth lone pair and oxygen p-orbital was crucial for short range structure formation.It affected macroscopic electrical properties including dielectric permittivity and ferroelectric behaviors, as mentioned later.The different characteristics observed within the band between 200-450 cm −1 and 450-700 cm −1 are likely associated with the change in local structure due to the presence of different A-cation combinations 18 .This indicated that local Ti-O bond and TiO 6 octahedral vibration of the NB, KB, NL, and KL ceramics was affected by the composition of the A site, particularly mono-/trivalent cations.The characteristic of the cubic phase, which is identified by a low intensity peak at around 800 cm −1 , was observed in NB, KB, and KL ceramics.Nonetheless, the NL peak exhibits an increase in intensity, confirming the transition from cubic to tetragonal phase [20][21][22][23][24] .It should be noted that the overall Raman spectra of bismuth containing compounds (NB and KB) exhibited similar peak shape and modes to BNT ceramic while the NL and the KL revealed strong deviation from the BNT ceramic.

Linear shrinkage, density, and relative density measurement
As shown in Fig. 2a-c, the linear shrinkage, density, and relative density of the sintered NB, KB, NL, and KL ceramics were evaluated to determine the effect of various mono-and tri-valent A cation substitutions on densification.The KB ceramic showed the largest linear shrinkage (16.1%).Meanwhile, the linear shrinkages of NL (4) and NB ceramics were lower with values of 14.91% and 13.92%, respectively.The KL ceramic exhibited the lowest linear shrinkage at 11.61%.Under same conditions, the NB and KB ceramics attained high densities of 5.28 g/ cm 3 and 5.26 g/cm 3 , respectively.The NB ceramic showed a slightly higher relative density than KB ceramic with values of 97.72% and 96.61%, respectively.The relative density value of 95% or more was expected for dense ceramics.This lower value was an indication of greater porosity 25 .In comparison, the NL and KL ceramics persisted lower densities with values of 4.77 g/cm 3 and 4.07 g/cm 3 , respectively.Both NL and KL ceramics were highly porous with the relative densities of 90.43% and 73.33%, respectively.This demonstrates that the ceramic was dense with compositions containing the Bi 3+ ion, as compared to the La 3+ ion.Since La 2 O 3 exhibited a melting point (2315 °C) which is significantly higher than that of Bi 2 O 3 (817 °C).Compositions containing lower www.nature.com/scientificreports/melting temperature components were able to densify at lower temperatures [26][27][28] .It should be noted that higher sintering temperature, i.e., 1275 °C, for the KL and NL samples were carried out.However, the obtain ceramics, which showed similar densities and shrinkage, exhibited partially melt surface and warping.Therefore, it was not suitable for further investigation of these samples.

Scanning electron microscope (SEM) microstructure analysis
Microstructure of NB, KB, NL, and KL ceramics sintered at 1250 °C were evaluated.As shown in Fig. 3a-d, the grain structure of NB ceramic was highly dense.The grain size as measured was 3.10 μm.The KB ceramic exhibited a dense and smaller grain size of 1.34 μm.The NL ceramic featured cubic-shaped grains.It exhibited grain size of 1.39 μm and showed visible porosity.The KL ceramic exhibited agglomerated particles with large pore sizes.All ceramics corresponded to the relative density values as illustrated in Fig. 2c.The grains in the NB and KB ceramics were larger and denser than those in the NL and KL ceramics.This could be due to the influence of Bi 3+ , which rapid diffusion coefficient of Bi 3+ ion allowed diffusion to occur more efficiently, thus grain growth was expected 26 .Meanwhile, the addition of La 3+ resulted in a reduced sintering diffusion rate, resulting in finer grains.It had been observed that the La 3+ ion could reduce the resultant grain size at low doping concentrations (≤ 0.2 at.%) by inhibiting grains growth [29][30][31] .Furthermore, the estimated elemental compositions of the ceramics were determined using XRF analysis.Table 2 lists the weight and atomic percentages of the elements, excluding oxygen, found in NB, KB, NL, and KL ceramics.The ideal atomic percentages should consist of 4 at% of each element on the A-site and 20 at% of Ti for compositions containing five equimolar elements on the A-site.The results indicate that atomic percentages of the constituent elements were estimated to be relatively close to their expected values from the batch calculations.It should be noted that the errors obtained were found to be less than 1.00% for Ti and Sr, and between 1.00 and 2.00% for other elements except light element.The higher error values of Na and K were from weak fluorescence signals.

Dielectric properties
The temperature dependence of the dielectric response was measured over the frequency range of 1 to 100 kHz and over the temperature range from room temperature to 350 °C.The relative permittivity and dielectric loss of the NB ceramics were shown in Fig. 4a.The maximum relative permittivity and dielectric loss were seen near Table 1.The details of refinement parameters.www.nature.com/scientificreports/room temperature and became broad towards higher temperatures.The relative permittivity tended to decrease with increased temperature, whereas the dielectric loss was relatively low and increased with increased temperature.It could be noted that the NB ceramic exhibited relaxor ferroelectric behavior from its dielectric response.

Composition
The relative permittivity and dielectric losses of the NB and KB ceramics were relatively similar values, 3000 and 0.15, respectively at room temperature.The temperature dependence of relative permittivity and dielectric loss of KB ceramic was displayed in Fig. 4b.A high relative permittivity peak gradually transformed to a broad peak over the wide temperature range.The similar frequency phenomenon appeared with the NB ceramic.This was an indicative of a relaxor ferroelectric behavior.Figure 4c revealed that the relative permittivity and dielectric loss of the NL ceramic decreased with increasing frequencies.The relative permittivity was about 400 to 200 from room temperature to 350 °C with dielectric loss of less than 0.2 over a wide range temperature.The drastic reduction in the relative permittivity as a function of temperature caused by the increase of degree of disorder of dipoles 32 .However, at high temperatures, from 250 °C to 350 °C, the dielectric loss increased gradually at higher temperature, suggesting the onset of conduction.This phenomenon revealed that relaxor ferroelectric behavior was not observed for the NL ceramic.Figure 4d showed a strong frequency dependence of relative permittivity and dielectric loss of the KL ceramic at the temperature of below 100 °C.At higher temperature, the relative permittivity and dielectric loss of the KL ceramic with various frequencies remained stable.However, there was an increase in dielectric loss of the KL ceramic at high temperature due to conductive characteristics.
To quantitatively characterize the relaxor behavior of the ceramics, the modified Curie-Weiss law was used to calculate the degree of relaxation characteristic or the degree of diffuseness (γ) which could be expressed by the following equations: where γ represented the dielectric relaxation factor which lie between 1 and 2, and ε m represented the maximum dielectric constant.The value of γ equal to 1 represents a normal ferroelectric transition whereas a γ value equal to 2 represents a relaxor ferroelectric transition 33 .The relationship between ln(1/ε r − 1/ε m ) and ln(T − T m ) of the ceramics was displayed in Fig. 5a,b.The γ value for the NB and KB ceramics was calculated to be 1.70 and ( 5) www.nature.com/scientificreports/1.75, respectively.These values were in a similar range of (1.46-1.85)as found in other titanate-based relaxor ferroelectrics 7,9,20,34 .The KB ceramic was found to have the maximum γ and exhibited a high dielectric permittivity of over 3000 at 1 kHz over a wide temperature range due to the presence of the relaxor state 35 .It should be noted that the dielectric characteristics of KL and NL did not show relaxor ferroelectric characteristics, i.e., dielectric peaks occurred at different frequencies.Therefore, the degree of diffuseness of the KL and NL sample was not analyzed.

Ferroelectric behavior and energy storage properties
As shown in Fig. 6a-d, the room-temperature P-E hysteresis loops under different applied electric fields at 1 Hz were demonstrated.It could be seen that the slim loops were observed in NB, KB, and NL ceramics.However, the rounded loops were exhibited in the KL ceramic, suggesting high loss behavior under high electric field.The NB, NL, and KL ceramics reached breakdown electric fields of around 80 kV/cm, whereas the KB ceramic could withstand only 50 kV/cm.The slim characteristics of the P-E loops is indicative of relaxor ferroelectric nature for the NB and KB ceramics.The NL ceramic showed linear dielectric behavior.The KL ceramic showed lossy dielectric behavior.The P-E loops measured at 50 kV/cm of ceramics were selected for comparison, and their corresponding electrical properties were presented in Fig. 7a-c.As for slim loops, the KB ceramic clearly displayed the greatest P MAX , followed by the NB and NL ceramics.All ceramics exhibited noticeable remnant polarization (P REM ) values  that were not drastically different.The highest P MAX of 16.9 μC/cm 2 with P REM of 1.23 μC/cm 2 were achieved in the KB ceramic.The NB ceramic exhibited P MAX of 13.1 μC/cm 2 with P REM of 0.93 μC/cm 2 .The NL ceramic showed P MAX of 1.83 μC/cm 2 with P REM of 0.84 μC/cm 2 .Coercive filed (E C ) values for all ceramics were in the range of 2.05 to 2.98 kV/cm for NB and KB ceramics.According to the concept of high entropy oxides, reduction of P REM resulted from random occupancy of multiple ions at the A-site, which disrupt polar ordering.Therefore, the NB and KB ceramics showed slim P-E loops with high P MAX and low P REM values, indicating relaxor behavior.The NL ceramic was a linear dielectric due to its linear response and lower dielectric permittivity.
The corresponding energy storage properties were calculated and plotted in Fig. 7d,e.Maximum recoverable energy density (W REC ) of 0.34 J/cm 3 was achieved in KB ceramic when applied to electric field of 50 kV/cm.The NB ceramic exhibited W REC of 0.27 (0.52 J/cm 3 at 80 kV/cm), whereas the NL ceramic showed W REC of 0.01 J/cm 3 .All ceramics exhibited similar energy loss density (W LOSS ) values.W LOSS for KB ceramic was 0.08 J/cm 3 , W LOSS for NB ceramic was 0.08 J/cm 3 (0.23 J/cm 3 at 80 kV/cm), and W LOSS for NL ceramic was 0.05 J/cm 3 .Consequently, the greatest energy efficiency (ŋ) was 81.51%, achieved by the KB ceramic.Lower value for NB ceramic was 78.65% (69.76% at 80 kV/cm), whereas those for NL ceramic was 24.99%.Higher energy storage density benefited from the relaxor behavior of KB and NB ceramics.Thus, it was more effective in storing energy due to the low energy barrier of nanodomains (PNRs).The energy storage density of the ceramics followed similar trends as other high entropy perovskite oxides based on titanate, such as (Na 0.2 Bi 0.2 Ba 0.2 Sr 0.2 Ca 0.2 )TiO 3 (1.02J/cm 3 at 145 kV/cm) 8 , (Bi 0.2 Na 0.2 Ca 0.2 Sr 0.2 Ba 0.2 )TiO 3 (1.32J/cm 3 at 110 kV/cm) 7 , (Bi 0.2 Na 0.2 Ba 0.2 Sr 0.2 Ca 0.2 )TiO 3 (1.3J/cm 3 at 150 kV/ cm) 9 , (Ba 0.2 Na 0.2 K 0.2 La 0.2 Bi 0.2 )TiO 3 (1.062J/cm 3 at 100 kV/cm) 36 , and (Bi 0.2 Na 0.2 K 0.2 Ba 0.2 Ca 0.2 )TiO 3 (0.684 J/ cm 3 at 129 kV/cm) 10 .Therefore, the KB ceramic was a promising candidate for using as energy storage devices.

Temperature-dependent ferroelectric behavior
Figure 8a,b showed the temperature dependence of P-E hysteresis loops of NB and KB ceramics at 50 kV/cm, respectively.At 30 °C (room temperature), the P-E loops of the NB ceramics displayed a relaxor ferroelectric behavior, and at temperatures between 30 °C and 90 °C, the P REM values tended to rise, which could be due to the increase of leakage current leading to increasing polarization current for relaxor ferroelectric 37 .When the temperature exceeded 100 °C, the P-E loops clearly exhibited a higher loss tendency.The P-E loops of KB ceramics showed the relaxor ferroelectric behavior at 30 °C.The P-E loops displayed a slim characteristic with decreased P REM , and remained relaxor ferroelectric as the temperature increased.This suggests diffuse phase transition behavior, which was typical of a relaxor ferroelectric material.The decrease in polarization exhibited a strong ergodicity 38 .Thus, the ferroelectric behavior of NB and KB ceramics was affected by the difference in monovalent cations and temperature.

Impedance and electric modulus measurement
The complex impedance and electric modulus measurements were conducted for the purpose of investigating the conduction mechanisms in these materials, with the goal of differentiating between bulk and grain boundary contributions to the relaxation phenomena.The temperature dependence of the complex impedance of NB and KB ceramics over the temperature range 400 °C-500 °C is shown in Fig. 9.At the measured temperatures, the data for the NB ceramic in Fig. 9a showed two semicircles and tended to decrease as the temperature increased.These features were attributed to the effects of grain (high frequency) and grain boundary (low frequency) contributions to the impedance.This could likely be due to the inhomogeneous distribution defects associated with volatile cations in the ceramic 39 .In contrast, the KB ceramic in Fig. 9b displayed a single semicircle and decreased as the temperature increased, indicating a decrease in the DC resistivity 40 .Detailed analysis on the trend of Z′′ and M′′ values for NB and KB ceramics as a function of frequency over the temperature range of 400 °C-500 °C was further analyzed to identify conduction mechanisms and relaxation process as shown in Fig. 10.The semi-log plots of normalized Z′′ and normalized M′′ as a function of frequency showed multiple peaks that relate to the transport and relaxation processes 40 .In Fig. 10a, the NB ceramic exhibited two peaks of normalized Z′′.The first peak (Z1) at lower frequencies was poorly resolved and largely disappeared as the temperature increased.A major peak (Z2) was apparent at all temperatures and clearly shifted to higher frequencies as the temperature increased.In the modulus data for NB at low temperatures, there is evidence of a small peak (M1) which was correlated to the main impedance peak (Z2).A much larger peak in modulus (M2) was apparent at all temperatures and was not correlated to the impedance peaks Z1 or Z2.This is indicative of an electrical heterogeneity in the microstructure which can most likely be attributed to differences in the stoichiometry of grain and grain boundary regions.In contrast, the impedance and modulus data for the KB ceramic showed close correlation with a large peak in impedance and modulus (Z3/M3) appearing at the same frequency at all temperatures.This is indicative of a more electrically homogeneous microstructure.
The variation of relaxation time (τ) as a function of temperature for NB and KB ceramics derived from the impedance data were shown in Fig. 11a,b, respectively.The τ value of both ceramics followed the Arrhenius relationship as presented in the Eq. ( 6) where τ o is pre-exponential factor, k B is Boltzman constant, T is the absolute temperature, and E a is activation energy.The activation energies obtained from this fitting are listed in Table 3.These data show that the NB ceramic had the higher values of E a in grain boundary than that of the bulk regions.For KB ceramic, the E a values obtained from Z″ and M″ were nearly similar which again is indicative of an electrically homogeneous structure.However, the activation energies in grain boundary and grain of KB ceramic were higher than in NB ceramic.As the Z′′ and M′′ peaks represented dielectric relaxation (localized conduction) and long-range conduction (non-localized conduction), respectively 41 .In other words, for NB ceramics, the separation of peak frequencies between Z′′ and M′′ indicated that the domination of relaxation process by localized conduction of charge carriers.On the other hand, the non-mismatch of peak frequencies between Z′′ and M′′ of KB ceramics suggested the non-localized conduction process of charge carriers.Overall, the selected titanate perovskites based on equimolar five-component A-cations could be prepared by conventional solid state reaction method.It was clear that monovalent cations (Na + or K + ) show slight effects on average and local structure formation as well as dielectric properties.More importantly, the compound consisting distinct trivalent cations, i.e., Bi 3+ or La 3+ , resulted in different local structure formation, dielectric properties, and ferroelectric behavior.This is likely due to the presence of the Bi 3+ ion with lone-pair electronic structure which drives structural distortions 42 .Thus, the polar structure of the nano scale (i.e.PNRs) could be maintained in the NB and KB ceramics.

Conclusion
High-entropy perovskite oxides composed of (A′ 0.2 A″ 0.2 Ba 0.2 Sr 0.2 Ca 0.2 )TiO 3 (A′ = Na + , K + , A″ = Bi 3+ , La 3+ ) were synthesized using a conventional solid-state reaction.The average grain sizes ranged between 1.34 µm to 3.10 µm and the (Na 0.2 Bi 0.2 Ba 0.2 Sr 0.2 Ca 0.2 )TiO 3 and (K 0.2 Bi 0.2 Ba 0.2 Sr 0.2 Ca 0.2 )TiO 3 ceramics exhibited high relative density of more than 95%.These compounds had the highest dielectric constant and the lowest dielectric loss of ~ 3000 and ~ 0.1 at 1 kHz, respectively.The (K 0.2 Bi 0.2 Ba 0.2 Sr 0.2 Ca 0.2 )TiO 3 ceramic exhibited the highest maximum polarization and lowest remnant polarization, corresponding to relaxor ferroelectric behavior.Additionally, it exhibited a recoverable energy density of 0.34 J/cm 3 and an energy storage efficiency of 81.51% at 50 kV/cm.The P-E loops of (Na 0.2 Bi 0.2 Ba 0.2 Sr 0.2 Ca 0.2 )TiO 3 ceramic displayed a larger loss when the temperature went above 100 °C.While the (K 0.2 Bi 0.2 Ba 0.2 Sr 0.2 Ca 0.2 )TiO 3 ceramic maintained its relaxor ferroelectricity and showed a slim characteristic with decreasing P REM .These findings suggest that the high-entropy titanate-based oxide containing K/Na and Bi at A-site can be synthesized into dense ceramics with a high dielectric permittivity and relaxor ferroelectric behavior.On the other hand, the compounds with La 3+ at A-site exhibit less desired properties due to the lack hybridization of 6s 2 of bismuth and oxygen p-orbital.Thus, this study offers another approach to develop highperformance relaxor ferroelectrics with the desired properties.

Figure 5 .
Figure 5. Relationship between ln(1/ε r − 1/ε m ) and ln(T − T m ) measured at 100 kHz of (a) NB and (b) KB ceramics.The red solid lines are the fit to the Modified Curie-Weiss law.

Figure 7 .
Figure 7. (a) P-E hysteresis loops, (b) relation between P MAX , P REM and composition, (c) relation between E C and composition, (d) relation between W REC , W LOSS and composition, and (e) relation between ŋ and composition of NB, KB, NL, and KL ceramics at 50 kV/cm.

Figure 8 .
Figure 8. Temperature dependence on P-E hysteresis loops of (a) NB and (b) KB ceramics measured under an electric field of 50 kV/cm and a frequency of 1 Hz.

Figure 9 .
Figure 9. Complex impedance plot between the imaginary part (Z″) and the real part (Z′) of (a) NB and (b) KB ceramics in the temperature range of 400 °C-500 °C.

( 6 )Figure 10 .
Figure 10.Frequency dependence of the imaginary part of impedance (Z″) and electrical modulus of (a) NB and (b) KB ceramics in the temperature range of 400 °C-500 °C.

Figure 11 .
Figure 11.Arrhenius plot of relaxation time (τ) for (a) NB and (b) KB ceramics in the temperature range of 400 °C-500 °C.

Table 2 .
XRF analysis of elemental compositions in the NB, KB, NL, and KL ceramics.

Table 3 .
Parameters obtained from fitting relaxation time (τ) of imaginary part of impedance (Z″) and Imaginary part of electric Modulus (M″).